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The Duffing Oscillator 


A Duffing equation describes any oscillator featuring a cubic stiffness term, i.e. a 
nonlinear elasiticity, such as 


#+6¢+ ax + Bx? =0. 


It is named after Georg Wilhelm Christian Caspar Duffing (1861 — 1944), a German 
engineer and inventor who wrote a seminal work on forced oscillations in 1918 (cf. 
[DUFFING 1918]). 

This oscillator exhibits a nice chaotic behaviour when driven by an external forcing 
function i.e. in the form of 


£+ 62+ ax + Bx® = ycos(uyt). 


Thorough theoretical treatments of this oscillator may be found in [KORSCH et al. 2008], 
[HOHLER 1988], and [CHANG 2017]. Apart from its interesting mathematical prop- 
erties it is a nice dynamic system which can be implemented on an analog computer 
without complicated scaling and which invites to playing with parameters to achieve 
all kinds of nice phase space plots. 

Figure 1 shows the implementation of the forcing function f(t) = ycos(wot). To 
achieve a minimum of harmonic distortion, the two amplitude limiting Zener diodes are 
connected to an integrator input with weight 1. ¢ introduces a tiny positive feedback 
signal to avoid a decreasing amplitude, while 7 controls the amplitude of the output 
signal. 

Figure 2 shows the straight-forward implementation of the DUFFING oscillator. 6 
was scaled down by a factor of a to allow for values of 3 up to 10. 
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Figure 1: Forcing function f(t) 
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Figure 2: DUFFING oscillator 





Dr. BERND ULMANN, Issue #29, 02-APR-2021 


ANALOG . ; 
PARADIGM * Analog Computer Applications 








A good parameter set to start from is shown in the following: 


a=1 
p=) 
y=1 
6 = 0.02 
wo =1 


€ = as small as possible 


Varying the parameters manually shows the very rich behaviour of this forced oscil- 
lator which even shows chaotic characteristics. The effect of wo is quite distinct. Two 
typical phase space plots generated by plotting x against —x are shown in figure 3. 
To get a wider range for y the output f(t) may also be connected to an input with 
weight 10 on the DUFFING oscillator sub-circuit. 
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Figure 3: Two phase space plots showing the chaotic behaviour of the DUFFING oscillator 
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